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Abstract 

We report on measurements of the angular dependence of the irreversibility 

temperature T, rr (6) in YBCI2CU3O7-8 thin films, defined by the onset of 

a third harmonic signal and measured by a miniature Hall probe. From the 

functional form of Tj rr (6) we conclude that the origin of the irreversibility line 

in unirradiated films is a dynamic crossover from an unpinned to a pinned 

vortex liquid. In irradiated films the irreversibility temperature is determined 

by the trapping angle. 
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I. INTRODUCTION 



The origin of the irreversibility line (IRL) in the field-temperature (H — T) phase diagram 
of high-temperature superconductors (HTS) is intriguing and still a widely discussed topic 
[.'Q-0|. Experimentally, this line is defined as the border line at which the magnetic response of 
the sample changes from irreversible to reversible. In HTS, large fluctuations and relatively 
weak pinning lead to a rich H — T phase diagram with a variety of dynamic and static 
transitions which can be responsible for the appearance of magnetic reversibility P-|5|JT(]|-P^| . 
Thus, a thorough experimental investigation of the IRL is important for the understanding of 
the vortex-lattice behavior in superconductors in general and of the mechanisms responsible 
for the onset of irreversible magnetic response, in particular. 

Several models, like thermally activated depinning P,|^J^, |l3| , |il|1 , vortex- lattice melting 
TSHHl and a transition from vortex glass to vortex fluid p3[- p5fl , were proposed to identify 



the origin of the IRL in HTS. Also, attention was given to the possibility of pinning in the 
vortex-liquid phase P, [10| , |ll| , p6| and to different dissipation mechanisms above the melting 
line P Jl3|j27| -p9| . Irreversibility due to geometrical |7||| or surface barriers [2^] have also been 
proposed, but these mechanisms are less probable in thin Y Ba2Cu^,0-js films with strong 
pinning. The irreversibility line may be affected by sample-dependent properties such as the 
nature and density of pinning centers and by intrinsic or extrinsic anisotropy. For example, 
in superconductors with columnar defects, the irreversibility line may either be identified 
with the Bose-glass transition [^|,p9H37|], or related to the concept of a trapping angle fl38f . 
The configuration of the columnar defects is also very important, since it affects possibility 
for different types of depinning mechanism. A splayed configuration, for example, inhibits 
creep from columnar defects |39| , |40[| . Similarly "crossed" defects (i.e. defects at angles ±9) 
were shown to act collectively, i.e., they introduce unidirectional anisotropy such that the 
current density reaches its maximum for magnetic field directed in a mid angle between 
defects 1I],[42|. 

Experimentally, the situation is even more complex, since different techniques (magne- 
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tization loops, field-cool vs. zero-field-cool DC magnetization, peak in the imaginary part 



of the first harmonic etc.) yield different IRLs To a great extent, the reliability of 

the determination of IRL depends on the criterion for the onset of the irreversibility. We 
determine the irreversibility temperature at given DC field by the onset of third harmonic 
in the ac response, which, we believe, is one of the most reliable methods for contactless 



determination of the IRL [44|. In most experiments Tj rr is measured as a function of the 
external field H. This information is insufficient to distinguish between different models for 
the origin of the irreversibility. Additional information, like the frequency dependence of 
the IRL |2^j45[1 or its angular variation [jl8 , 2~0|, 21], 



is needed. 



In this paper we report on a study of the angular dependence of the irreversibility tem- 
perature T irr (8) in thin Y Ba^Cu^O-j^^ (YBCO) films before and after irradiation with Pb 



ions. 



II. EXPERIMENTAL 

The 1500 A YBCO films were "sandwiched" between SrTiOs layers P7| . First, a 500 A 
layer of SrTiO^ was deposited on a MgO substrate. Then, the YBCO film was laser 
ablated on top of the SrTiO^ and finally, the YBCO was covered by a protective 300 A 
layer of SrTiO^. All three samples have the same lateral dimensions of 100 x 500 /im 2 . 
One film, denoted as REF, was used as a reference sample. The other two, UIR and CIR, 
were irradiated at GANIL with 2 x 10 11 ions/cm 2 5.8 GeV Pb ions along the c— axis and 
along 9 = ±45°, respectively. (UIR and CIR stand for "uniform irradiation" and "crossed 
irradiation", respectively). The superconducting transition temperatures, measured by a 
Quantum Design SQUID susceptometer and defined as the onset of the Meissner expulsion 
in a DC field of 5 G, are T c « 89 K for the samples REF and UIR and 88 K for CIR. 

For the ac measurements reported below we used a miniature 80 x 80 \im 2 InSb Hall- 
probe, which was positioned in the center of the sample. The 1 G ac magnetic field, always 
parallel to the c— axis, was induced by a small coil surrounding the sample. An external dc 
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magnetic field, up to H a = 1.5 Tesla, could be applied at any direction 9 with respect to 
the c— axis. In our experiments DC magnetic field was always turned on at a fixed angle 
at T > T c and then the ac response was recorded during sample cooling. The irreversibility 
temperature, Tj rr (9) is defined as the onset of the third harmonic signal in the ac response 
measured by the Hall probe This procedure was repeated for various dc fields and at 
various angles 9 of the field with respect to the c— axis. 



III. RESULTS 

Fig. m presents measurements of V3, the third harmonic in the ac response, versus tem- 
perature T, during field-cooling at 1 Tesla for the sample REF at various angles between 
and 90°. Apparently, as the angle 9 increases the whole V3 curve shifts to higher tempera- 
tures and becomes narrower. The onset of irreversibility, T irr (9), is defined by the criterion 

yonset = q q 5 j n the unitg Q f p ig g 

Fig. |2] exhibits typical Ti rr (9) data for the unirradiated sample REF, measured at two 
values of the external field: 0.5 Tesla and 1 Tesla. Both curves exhibit a shallow minimum 
around 9 = and they reach their maximum values for H along the ab plane, at angles 
9 = ±90°. We also measured the frequency dependence of T irr for the same values of 
magnetic field. As shown in Fig. the slope dT irr /dln (/) is larger for larger field. 

The sample irradiated along the c— axis exhibits additional feature - a peak around 9 = 0. 
This is clearly shown in Fig. |] where we compare Ti rr (9) at H = 1 Tesla for the samples 
REF and UIR. As discussed below, this peak is a signature of the unidirectional magnetic 
anisotropy induced by the columnar defects. Intuitively, one would therefore expect two 
peaks, along 9 = ±45°, for the third sample, CIR, crossed-irradiated at 9 = ±45°. Instead, 
we find one strong peak around 9 = 0, similar to that found in BSCCO crystals pi2| . This is 
demonstrated in Fig. |5| where we compare T{ Tr (9) at H = 0.5 Tesla for this sample (CIR) 
and for the unirradiated sample (REF). We argue below that the peak around 9 = 0° is a 
result of a collective action of the crossed columnar defects, and its origin is the same as 
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that for unidirectional enhancement of critical current density observed in BSCCO crystals 



IV. ANALYSIS 

The "true" irreversibility temperature To is defined as a temperature above which the 
pinning vanishes. Such disappearance of pinning can be of static (true phase transition), as 
well as of dynamic origin (gradual freezing, pinning in liquid). In practice, one determines 
the irreversibility temperature T irr (A) as the temperature above which the critical current 
density is less than some threshold value A. Therefore, by definition, T = hm (T irr (A)). 
The apparent current depends on temperature T, magnetic field B and the frequency / of the 
exciting field which defines a characteristic time-scale 1/f for the experiment. By solving 
the equation j (T, B, f) = A with respect to T one finds the experimental irreversibility 
temperature T irr for constant B and /. In the following we argue that in our experiments 
the measured T irr is a good approximation of T . In order to estimate T irr we employ a 
general form for the apparent current density in the vicinity of the irreversibility line (IRL) 
Mi: 



;i - T/T ) Q / / 
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j (T, B, f) oc j c (0) v , ' MM (!) 



(B/B f V/c 



o , 



where the parameters B Q and fo are temperature independent (Eq. [3] is thus valid only in a 
narrow temperature interval near the IRL and for fields larger than H c i). From Eq. p] we 
get: 

T,r=n(B)\l-(-$-(i-)U^Y) I (2) 



^(0) \BJ \f 

Inserting reasonable numerical estimates: j c (0) ~ 10 7 A/ cm 2 , A ~ 100 A/ cm? for our 
experimental resolution, B ~ 10 3 G, B ~ 10 4 G, (3 ~ 1 @, 7 ~ 1 @, / ~ 10 2 Hz, 
and f ~ 10 7 Hz 0, we get from Eq. |: T irr = T (B) (l - O.OOS 1 /"). Thus, with 0.5% 
accuracy we may say that T, rr , the measured onset of the third harmonic component in the 
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ac response, marks some "true" irreversibility crossover line To (B). The nature of this line 
T (B) is our main interest, since, as discussed in the introduction, it is directly related to 
the pinning properties of vortex lattice in type-II superconductors at high temperatures. 



A. Unirradiated YBCO film 

We turn now to consider the effect of the intrinsic anisotropy on T rr (8). Following the 
anisotropic scaling approach proposed by Blatter et al. [B^,B0], we replace T by eT and B by 



B e ff = EgB, where eq = ycos 2 (6) + e 2 sin (6) and e ~ 1/7 is the anisotropy parameter for 
YBCO. It should be emphasized that we can use this scheme only in the case of intrinsic 



anisotropy e = ym a b/m c , where m c and m a fe denote the effective masses of the electron 
along the c— axis and in the ab— plane, respectively. In the case of some extrinsic magnetic 
anisotropy, (columnar defects or twin planes), the critical current depends on the angle not 
only via the effective magnetic field -B e //> but also because of this extrinsic anisotropy. 

As we have already indicated in the Introduction, there are several possible origins for 
a crossover from irreversible to reversible magnetic behavior in unirradiated samples. We 
exclude the vortex-glass to vortex fluid transition as a possible origin for the IRL, because this 



transition was shown to occur at temperatures lower than the onset of dissipation p4| , |26| , |5T 
The thermal depinning temperature increases with increase of field H T^ p oc \f~B and, 
therefore is excluded as well. Vortex-lattice melting transition is believed to be responsible 



for the appearance of reversibility [|T^,|T(||T^,|T!| . The explicit angular dependence of T m was 



derived by Blatter et al. [H,fE| using their scaling approach: 



T m W , 2 « ( W ^ « (l - |) ) (3) 

where $o is the flux quantum, £ is the coherence length, (3 m w 5.6 is a numerical factor, 
estimated in |J, cl — 0.1 is the Lindemann number, Gi = (T c /eH c2 (0) £ 3 (0)) 2 /2 is the 
Ginzburg number, and H c2 (0) is the linear extrapolation of the upper critical field from T c 
to zero. Solving Eq. § with respect to T m we get: 
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Eq. (| predicts that the melting temperature decreases as B e ff increases. This is due to 
the fact that the inter- vortex distance Oq oc 1/B e ff decreases faster than the character- 



istic amplitude of fluctuations (u 2 (B e ff,T m )) th oc l/y5 e //. Therefore, the condition for 
the vortex-lattice melting (u 2 (B e ff,T m )) th ~ c|oq implies larger melting temperatures for 
smaller effective fields, i. e., for larger angles. In agreement with this prediction, the ex- 
perimental data of Fig. show that Tj rr increases with the angle, i.e., decreases with 
B e ff. The solid lines in Fig. ^| are fits to equation |j. From this fit we get C ~ 0.0005. 



However, a reasonable estimate of C ~ W (j3 m Gi/ 'c\H C 2 (0)) yields C ~ 0.01, where we take 
H c2 (0) = 5 ■ 10 6 G, c L = 0.1, Gi = 0.01, and (3 m = 5.6 §. Also, Yeh et al. showed that the 
onset of irreversibility occurs above the melting temperature (Ref. ||29|| , Fig. 4). In addition, 
the important effect of the frequency (see Fig. |^) is not included in Eq. f|. 

We discuss now another possibility for the onset of the irreversibility, namely, pinning in 
the vortex liquid (for a discussion see Ch. VI in Blatter et al. || and references therein). 
Any fluctuation in the vortex structure in the liquid state has to be averaged over the 
characteristic time scale for pinning t P i n . In the absence of viscosity the only fluctuations 
in the liquid state are thermal fluctuations, which have a characteristic time t t h << t P i n . 
(As shown in || t pin /t t h oc jo/j c , where j is the depairing current). Thus, such a liquid is 
always unpinned. The situation is different for a liquid with finite viscosity. In this case 
there exists another type of excitations in the vortex structure, i.e. plastic deformations with 
a characteristic time scale t p i. The energy barrier, corresponding to plastic deformation is 
shown to be P,|T2||: 



U pl ~ 7 eeoao * 7 (|§) ' (T c - T) B^ 2 (5) 
where 7 is a coefficient of the order of unity. The corresponding characteristic time scale is: 

t p i ~ t th exp(U p i/T) (6) 



Thus, depending on the viscosity, t p i can be smaller or larger than t pin . In the latter case, 
after averaging over a time t pin , the vortex structure remains distorted and such a liquid 
shows irreversible magnetic behavior. Thus, on the time scale of t P i n the distorted vortex 
structure is pinned. The crossover between pinned and unpinned liquid occurs at tempera- 
ture Tfc where the characteristic relaxation time for pinning t pin (T) becomes comparable to 
that for plastic motion t p \ (T). Thus, using Eqs. || and |B] we obtain: 

T k = (7) 

1 ( AGi ^ 2 In (t pm /t th )VB 



1 7 \n c2 (o) 

Finally, using the anisotropic scaling we may rewrite Eq. [7| for f pin < f < fth as: 

T irr (9) = T k (9) = -J* = Tc (8) 

with f t h = 1/tth and f p i n = l/t pin . Note the apparent similarity with the expression for the 

i 

melting temperature, Eq. [|. The numerical estimate for A = ^ (jfj^j 2 ^(fth/f) gives: 
A ~ 10 _4 ln (fth/ ' f) /t- This is in agreement with the value found from the fit (solid line in 
Fig. |) for H c2 (0) = 5 ■ 10 6 G, Gi = 0.01, f th ~ 10 10 Hz, and 7 ~ 4. 

To further confirm that in our YBCO films the most probable physical mechanism for 
the onset of irreversibility is a dynamic crossover from unpinned to pinned vortex liquid we 
discuss now the frequency dependence of T irr . Equation [8] has a clear prediction for the fre- 
quency dependence of Ti rr . To see it directly we may simplify it by using the experimentally 
determined smallness of value of the fit parameter A ~ 0.0005, which allows to expand Eq. 
H (for not too large fields) as 

AGi 



' ' ^ :/ ' ' 1 1 ~ 7 [h^W)J ln {fth/f) v eB 1 1 9 ] 

which results in a linear dependence of T irr upon ln(/) and a slope S = dTirr/d ln(/) ~ 
Y ( j/^o) J 2 V £ eB = T c A\JeqB ln (f / fth)- Note that the slope is proportional to This 
is indeed confirmed by the experimental data, as is demonstrated by the solid lines in Fig. 
I3L From this fit we get S/y/B = 0.004 and we can independently verify the parameter A 
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appeared in Eq. § A = S/ [T cV ^B\n (f/fth)) = 0.0008, which is in an agreement with the 
value obtained above. 

We note that the approximated expression for the frequency dependence of T irr , Eq. |9|, 
is valid in the whole experimentally accessible range of magnetic field since Eq. [5] predicts a 
maximum in the slope S at -B max = (Ay/£e} ~ 400 Tesla for the experimental parameters. 
This value is, of course, beyond the experimental limits and, probably, even exceeds H c2 . 

Another support for the onset of the irreversibility in a vortex liquid is the ac field 
amplitude dependence of the IRL. In both, thermal- activated (TAFF) and pure flux-flow 
(FF) regimes the I-V curves are linear and the onset of the third harmonic is due to a 
change in the slope (from p FF to Ptaff)- In this case we expect the amplitude dependence 
for this onset. Contrary, at the melting transition the onset of irreversibility is sharp and 
is not expected to depend upon the amplitute of the ac field. In our experiments we find a 
pronounced amplitude dependence of the IRL, thus confirming the above scenario. 

B. Irradiated YBCO films 

For the irradiated films the situation is quite different. The models for T irr (8) in unirra- 
diated films cannot explain the experimental features exhibited in Figs. |] and|^, in particular 
the increase in Tj rr in the vicinity of 9 = 0. Such discrepancy can only be due to the angular 
anisotropy introduced by columnar defects, i.e., the angle-dependent pinning strength. It 
was shown, both theoretically and experimentally ||37|| , that for magnetic field oriented 
along the defects the irreversibility line is shifted upward with respect to the unirradiated 
system. Thus, our results in Fig. [| suggest that the measured T irr (6) is a superposition of 
the angular variation of Tj rr in unirradiated film (denoted in this section as T^ F ) and the 
anisotropic enhancement of the pinning strength due to irradiation. 

We can estimate the latter contribution by employing the concept of a "trapping angle" 
6 t , the angle between the external field and the defects at which vortices start to be partially 
trapped by columnar tracks. (For a schematic description, see Fig. 43 in Blatter et al. ||). 
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As we show in Appendix A, 



tan(e t )^j2e r /e l (10) 



where e r (T) is the trapping potential of a columnar defect and E\ is the line tension. In 
the experiment we cool down at a fixed 9, and the onset of irreversibility must occur when 
9 = 9 t (T), provided that the temperature is still larger than T^ FF (9). Otherwise, the onset 
occurs at T^ FF . This defines the condition for the irreversibility temperature Tj rr for angles 
9<9 c = 9 t (t£ ff (0 t j) ~ 50° in our case: 



tan (9) = tan (9 t (T irr )) w ^2e r /e, (11) 

At high-temperatures e r (T) oc exp (—T/T^j, where T r dv is the depinning energy ||. Thus, 
we can write for Tj rr : 

rr \T^ EF (9)-Dln(C\Un(9)\) 9 < 9 C 

J-irr{v) = < (12) 

( T*** (9) 9 > 9 C 

where D and C are constants. This expression is in an agreement with our results shown 
in Fig. |] (solid line). We note, however, some discrepancy in the vicinity of 9 = 0, where 
we find quite weak dependence of T irr on angle. We explain this deviation by considering 
the influence of relaxation, which, in the case of parallel defects depends on angle. The 
relaxation rate is maximal, when vortices are aligned along the defects and retains its nor- 
mal "background" value for perpendicular direction []52| . Vortex, captured by a defect, can 
nucleate a double kink which slides out resulting in a displacement of a vortex on a neigh- 
boring column. In our irradiated samples the defect lattice is very dense (the matching field 
Bfj, = 4 Tesla, i.e., distance between columns d ~ 220 A) and such double-kink nucleation 
an easy process. Thus, the irreversibility temperature should be shifted down around 9 = 
as compared to the "ideal", non- relaxed value, Eq. This explains the reduction in T irr 
in Fig. f§ 

We may now conclude that in irradiated films, for angles less than the critical angle 9 C , the 
irreversibility line is determined by the trapping angle 9 t . The Bose-glass transition can 
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probably only be found for small angles within the lock-in angle 9l < 10°. This conclusion 
is also indirectly confirmed in |kJ. 

As was pointed out in the introduction, crossed defects should hinder the relaxation due 
to forced entanglement of vortices. Thus, the irreversibility temperature is expected to be 
closer to that predicted by Eq. |T2|. Fig. |5| shows a good agreement of the experimental data 
with Eq. [12] (solid line). To explain why defects crossed at large angle act collectively and 
force unidirectional magnetic anisotropy, we follow here the approach outlined in and 
extend that description to account for arbitrary orientation of the external field with respect 



to the crossed columnar defects and to the c— axis. In Ref. |4lJ the authors consider possible 
motion of vortices in a "forest" of crossed defects for field oriented along the c— axis. In 
our case of a dense lattice we may exclude from consideration free kink sliding and consider 
only depinning from the intersections. We sketch in Fig. |6| the two limiting situations: (a) 
the external field is parallel to one subsystem of the columnar defects (9 = 45°) and (b) the 
external field is oriented along the c— axis, between crossed columns (9 = 0). In case (a), Fig. 
^|a, vortices can depin just by nucleation the single kinks which are sliding from intersection 
to intersection, or, by nucleation of super-kinks resulting in a kind of motion, similar to a 
variable-range hopping. This type of thermally assisted vortex depinning does not cost any 
additional energy on vortex bending. Another situation arises for field along the c— axis, Fig. 
§3. Now vortices can depin only via nucleation of multiple half-loops, which characteristic 
size depends upon current density. This results in additional barrier for vortex depinning, 
which even diverges at zero current ||. As a result, the relaxation rate is anisotropic, i.e, 
it is suppressed when the external field is oriented along the mid direction between the two 
subsystems of the crossed columnar defects. This is just opposite to a situation in uniformly 
irradiated samples. 
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V. CONCLUSIONS 



We presented angle-resolved measurements of the irreversibility temperature in unirradi- 
ated YBa 2 Cu 3 7 _s fil m and in two films with columnar defects, induced by 6 GeV Pb-ions 
irradiation, either parallel to the c— axis or 'crossed' in = ±45°. We find that in the unir- 
radiated film the transition from irreversible to reversible state occurs above the melting 
line and marks the crossover from pinned to unpinned vortex liquid. In irradiated 
films, within the critical angle 9 C ~ 50°, the irreversibility line is determined by the tem- 
perature dependent trapping angle. For larger angles T irr is determined by the intrinsic 
anisotropy via the effective field. The formulae for T irr (9) for both unirradiated and irradi- 
ated films are given. We also discuss the possible influence of anisotropic enhancement in 
relaxation rate which leads to a smearing of the expected cusp at 9 = in the T irr {9) curve 
in the uniformly irradiated film. Finally, we demonstrate the collective action of crossed 
columnar defects, which can lead to suppression of relaxation and enhancement of pinning 
strength along the mid direction. 
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APPENDIX A 



We describe here the derivation of our Eq. 10, which differs slightly from the analogous 
Eq. (9.173) of Blatter et al. ||. We derive it using exactly the same approach (and notions) 
as in H, but, in view of the experimental situation, avoid assumption of small angles, which 
allows Blatter et al. to approximate tan (8) ~ sin (8) ~ 8. In order to estimate the trapping 
angle one has to optimize the energy change due to the vortex trapping by columnar defects. 
This energy is written as [H]: 



e (r, 8) = ei 



r + <T + 7— r - r 



d \\ d 



- re r (1A) 



tan (8) J / sin (9) 

where r (8) is the length of the vortex segment, trapped by a defect; d is the distance 
between the columns; E\ is the line tension; and e r is the trapping potential of the defects. 
The variation of Eq. |1A| with respect to r at fixed a angle 8 defines the angular dependence 
of r (8). The trapping angle 8 t can be found by solving the equation r (8 t ) = 0. This results 



111 



Je r (2ei — e r ) 



which, at sufficiently small e r , can be approximated as 



tan 



^ + 0(e 3 / 2 ) (3A) 



Apparently, at very small angles we recover the original result of Ref. M. In the paper, 



for the sake of simplicity, we use expression Eq. |3A| instead of the full expression Eq. |2A . 
However, as noted above we cannot limit ourselves bto small angles and, generally speaking, 
the trapping angle may be quite large (8 t w 40° in our case). The error due to use of Eq. 



3A" can be etimated as follows: at 8 ~ 40° Eq. gives e r /ei ~ 0.24, whereas Eq. |3A| gives 



e T /ei ~ 0.35, which is suitable for our implication of Eq. |3A| , since we consider exponential 
decrease of e r . Also, as shown in in a system with anisotropy e, the trapping angle is 
enlarged by a factor of 1/e. 
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FIGURES 

FIG. 1. The third harmonic signal V3 versus temperature during field-cooling at 1 Tesla for 
sample REF at 6 = 0, 10°, 30°, 40°, 60°, 80°, 90°. 

FIG. 2. The irreversibility temperature in the unirradiated sample REF at two values of the 
external field: H = 0.5 and 1 Tesla. The solid lines are fits to Eq. H 

FIG. 3. The frequency dependence of Ti rr in the unirradiated sample REF at two values of the 
external field: H = 0.5 and 1 Tesla. The solid lines are fits to Eq. |. 

FIG. 4. The irreversibility temperature for two samples: REF (unirradiated - open circles) and 
UIR (irradiated along the c— axis sample, - filled circles) at H = 1 Tesla. Solid lines are fits to 
Eq. |8] and Eq. [l^, respectively. 

FIG. 5. The irreversibility temperature for two samples: REF (unirradiated - open circles) and 
CIR (irradiated along 9 = ±45° - filled circles) at H = 0.5 Tesla. Solid lines are fits to Eq. || and 
Eq. [l2|, respectively. 

FIG. 6. Schematic description of a possible depinning modes of a vortex line in the case of 
crossed columnar defects; - (a) magnetic field is directed along 9 = 45°; (b) magnetic field is along 
9 = 0. 
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